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Abstract 

The most general nonuniform reaction-diffusion models on a one-dimensional 
lattice with boundaries, for which the time evolution equations of corre- 
lation functions are closed, are considered. A transfer matrix method is 
used to find the static solution. It is seen that this transfer matrix can be 
obtained in a closed form, if the reaction rates satisfy certain conditions. 
We call such models superautonomous. Possible static phase transitions 
of such models are investigated. At the end, as an example of superau- 
tonomous models, a nonuniform voter model is introduced, and solved 
explicitly. 
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1 Introduction 



Most of the investigations on reaction-diffusion models are devoted to uniform 
models, where interaction rates are site-independent. Among the simplest gen- 
eralizations beyond a completely uniform system is a lattice with alternating 
rates. In [TJ, relaxation in the kinetic Ising model on an alternating isotopic 
chain has been discussed. In [2HI], the steady state configurational probabilities 
of an Ising spin chain driven out of equilibrium by a coupling to two heat baths 
has been investigated. An example is a one-dimensional Ising model on a ring, 
in which the evolution is according to a generalization of Glauber rates, such 
that spins at even (odd) lattice sites experience a temperature T e (T a ). In this 
model the detailed balance is violated. The response function to an infinitesimal 
magnetic held for the Ising-Glauber model with arbitrary exchange couplings 
has been studied in [5]. Other generalizations of the Glauber model consist of, 
for example, alternating-isotopic chains and alternating-bound chains ( [3] for 
example). In a recent article [7], we studied the expectation values of spins in 
an Ising model with nonuniform coupling constants. A transfer matrix method 
was used to study the steady state behavior of the system in the thermodynamic 
limit. Different (static) phases of this system were studied, and a closed form 
was obtained for this transfer matrix. 

In [8] a ten-parameter family of one-species reaction-diffusion processes with 
nearest-neighbor interaction was introduced, for which the evolution equation 
of n-point functions contains only n- or less- point functions, the so called au- 
tonomous models. The average particle-number in each site was obtained ex- 
actly for these models. In [9j[10] , this was generalized to multi-species systems 
and more-than- two-site interactions. In |11IU3] . the phase structure of some 
classes of single or multiple-species reaction-diffusion systems was investigated. 
These investigations were based on the one-point functions of the systems. 

In the present paper the most general nonuniform exclusion nearest-neighbor 
reaction-diffusion models on a one-dimensional lattice with boundaries are stud- 
ied, for which the evolution equations of the one-point functions are closed, and 
the transfer matrix has a closed form. It is shown that there is a possible 
phase transition in such models, which corresponds to a reduction of the role 
of boundary conditions on time-independent profile of the expectation value of 
the number operators. The scheme of the paper is as follows. In section 2, the 
models are introduced, and the evolution equation for the expectation values 
of rij (the number operators at the site i) is obtained. Also conditions are ob- 
tained so that the evolution of the expectation values of 7ij is closed. In section 
3, the equation governing the static solution for the expectation values of n^s 
is obtained, and a transfer matrix method is introduced to obtain the static 
solution and investigate different (static) phases of the system. It is also seen 
that to write a closed form for the transfer matrix, further conditions on the 
reaction rates are to be satisfied. We call models satisfying these conditions su- 
perautonomous models. In section 4, as an example, a nonuniform voter model 
is investigated in more detail. Section 5 is devoted to the concluding remarks. 
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2 Exclusion nearest-neighbor reaction-diffusion 
models with nonuniform reaction rates 



Consider a one-dimensional lattice with (L + 1) sites, numbered from to L. 
Each site is either empty (denoted by the vector eo) or occupied with one particle 
(denoted by the vector e — 1). The evolution of the system is said to be governed 
by nearest-neighbor interactions, if the evolution of each site depends on only 
that site and its nearest neighbors (sites directly related to it through a link). 
The evolution of such a system is governed by a Hamiltonian % of the form, 

n = n' +^2n a ^+n' L , (i) 

where H a corresponds to the link a: 

H a = l®( Q -"> ® H a ® (2) 

and 

M=4 (3) 

The link a links the sites (a — [/,) and (a + fj,), so that a ± \i are integers, and a 
runs from ji up to (L — /i). Throughout this paper, sites are denoted by Latin 
letters which represent integers, while links are denoted by Greek letters which 
represent integers plus one half (/i), so that the link a joins the sites (a — fj,) and 
(a + fi), while the site i joins the links (i — fj,) and + H a is a linear operator 
acting on a four dimensional space (the configuration space corresponding to the 
sites (a — fi) and (a + fi)) with a basis {e o, eoi, eio, en}. Also, 

U' = H' Q ®l® L , 

U' L = 1® L ®H' L , (4) 

where H'O and H' L are linear operators acting on two dimensional spaces (the 
configuration spaces corresponding to the sites and L, respectively) with bases 
{e ,ei}. The nondiagonal components of H a , Hq, and H' L are reaction rates. 
Denoting a full site by • and an empty site by o, the possible reactions for the 
boundary sites or L are 

o — > •, with the rate {H^^) 1 ^, 

• -> o, with the rate (H' QL )°i, (5) 
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while those for the link a are 
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(6) 

It is seen that these rates are in general different for different links (also the rates 
on the boundary sites are in general different). The system is called uniform if 
the rates are the same for all links, and nonuniform if it is not the case. 
The number operator in the site i is denoted by nf. 



1®' <8> n ® l®^-*) 



(7) 



where n is an operator acting on a two dimensional space with the basis {e , e{\. 
The matrix form of n in this basis is 

n a b = StSl (8) 

The evolution equation for the expectation value of an observable Q is 

d 



where 



(9) 
(10) 

the vector >J is the (2 i+1 dimensional) probability vector describing the system 
and S is the covector 



dt (Q) = (QH), 
(Q) = SQV, 



and 



S := s 

So = 1- 



(11) 

(12) 



The system is called autonomous, if the Hamiltonian is so that the evolution 
of the expectation values of rii is closed in terms of the expectation values of rij 's. 
In the evolution equation for the expectation value of rij, the expectation values 
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of rij_i, rii, rii+i, (rii—ini), and (r^rii+i) occur. It is seen that the criterion 
that the coefficients of the last two vanish, is 

s a [H^r® r] al =0, 

s a [H i+ll r® r] la =0, (13) 

respectively, where 

r Q = -(52 + <^. (14) 

Equation ([13")) should hold for all z's, in order that the system be autonomous. 
So one can rewrite it like 

s a [H a r® r} al = 0, 

s a [H a r® r] la = 0. (15) 

It is seen that this condition is the same as the corresponding condition for 
uniform lattices [HUH], written for each link separately. Provided that this 
condition holds, one arrives at 

^("i) = Vi-n ("i+l) + («i-/i+I / t+A t ) + < i < L, 

(16) 



where 
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(17) 

For the boundary sites (the sites and L), one has 

±(no) = 9, <m> + - K)°i + ^] (no) + + (18) 

d_ 

dt' 

(19) 



:(n L ) = ( ni _i) + + (i^l) 1 ! - (^)°i] (n L > + + (H' L )\} 



3 The static solution 

For the static solution ((n) s t), the left hand side of p^|) vanishes and one obtains 
(n i+ i) st = (ni-x) st - (m)st -. (20) 

Vi+H "i+ti 
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Denoting that part of this solution which satisfies the homogeneous equation by 
(n)£° ro , it is seen that 



<m>rs 



/„, \hom Vi-fi i m \hom + V i+[J. i \hom 

\ n i+l/st = — a \"i-i/st a 

which can be written as the following matrix form 



horn 



rr = d, x 



where 



and 



Dt := 



1 



(21) 

(22) 
(23) 

(24) 



Using these, one arrives at 



Xhom y\ yhom 



(25) 
(26) 



where 

D a /3 '■= D a -n Da-fi-i ■ ■ ■ Dp +fl . 

To solve (|20)). one can use a tranfer matrix (Green's function) method. Consider 
the equation 



Gi+lj — . I,; 



Defining 



y — 



Ga-nj 



it is seen that the solution for Y is 



DapYpj, a < j 
DapYpj, a > j 



with the condition that (|2"T|) holds for i = j. This condition is 



which reads 



(27) 
(28) 

(29) 

(30) 
(31) 
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A particular solution for Y is obtained if one sets Y equal to zero. In this case, 



Y a j — * a j Daj+ft 

where O is the step function: 



(32) 



9 ai := I ' a< J.- (33) 

1, OL > J 



Using this, the general solution to (f2T)|) can be written as 



X a — D a p Xp orVL — D a j 



j<a 7 3+» 



(34) 



As it was the case in [7], the steady state profile near the end-site is 
determined by the eigenvalues of the matrix D a ^ , where a is some site far from 
the ends. One has 

X a = X® i a , 

A^ = X^ t a , (35) 

where f a is the eigenvector of D a ^ corresponding to the eigenvalue A°, and 
X°'s and I°'s are the coefficients of expansions of X a and X^ in terms of the 
eigenvectors. If there was no nonhomogeneous part in the equation (1161) . then 
the discussion would be exactly similar to that of 7 : A^ vanishes if A a tends 
to infinity (in the thermodynamic limit). Otherwise, X° is generally nonzero 
and determined by the boundary conditions. If the nonhomogeneous part does 
not vanish, then the second term in Q34[) generally contains a large multiple of 
f a if A a tends to infinity. This large part is to be cancelled by a large multiple 
of f a coming from the first term in (l34l) . So in this case X° does not vanish 
but tends to a fixed value independent of boundary conditions. It is seen that 
although equations (fl8|) and (|T9|) serve as boundary conditions to obtain say A M , 
the above general argument is independent of these conditions. The essence of 
the above argument is the following. In general, X a (the components of which 
are noting but expectations of number operators) is a linear combination of 
two vectors (f a 's) plus a nonhomogeneous part. The nonhomogeneous part is 
determined from the bulk reactions, and the two unknown coefficients f a 's are 
to be determined from the the boundary conditions resulted from (|18p and 
(Tl9|) (in the static case that the left hand sides vanish). There are, however, 
regions of the parameter space where in the thermodynamic limit one of the 
coefficients of f Q 's (or possibly both) are determined from the bulk reactions 
only. This essentially means that the effect of boundaries on the behavior of 
the system is reduced. In a transport system, for example, it is expected that 
the time-independent profile of the moving bodies' density depend on both the 
bulk reactions (speed, overtaking, etc) and the boundary reactions (injection 
and extraction rates) . But there could be cases where these boundary terms are 
unimportant, or less important. 

The situation (in the thermodynamic limit) can be summarized as follows. 
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i The eigenvalue A a tends to infinity. In this case 1° tends to a fixed value, 
independent of the boundary conditions. 

ii The eigenvalue A a tends to zero or a finite number. In this case X° is 
determined by the boundary condiotions. 

Obviously, similar cases occur at the other boundary site. It is seen that this 
behavior at one of the boundaries is independent of the analog behavior at the 
other boundary. 

This is the static phase transition of the system, some discontinuous behavior 
of the expectation value of the number operator near the boundaries (note that 
the components of X a are nothing but the expectations of number operators). 

If one can write Di as 



(36) 



where A* is diagonal, and E Q depends on only the parameters corresponding 
to the link a, then it is easy to find the solution to (j2"2")l . The system is called 
superautonomous, if this is the case. Putting 
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where 

s~a . a a d a b a c Q 
Using (|39f and (|40|) . one obtains A, and Be 

Ci 



A, 



Bi = 



^l — fl 

di—fj. 



(37) 
(38) 

(39) 
(40) 
(41) 
(42) 

(43) 

(44) 
(45) 
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Using these in (|41]h it is seen that 

— j — = 7—^, (46) 

which can be solved as 



a a b a = c/)8 a <, a , (47) 
c a d a = <i>T] a <; a , (48) 

where is a constant. Finally, using (J44J), (J45J), (|37)l . and (08]) in (j42|), on obtains 



1 = —q> Hi-fi 



$Vi+^ (49) 



the solution to which is 



a a d a = (-<pK a + ip)<; a , (50) 

b a c a = {-(j>v a -ip)q a , (51) 

where ip is another constant. <; is not an independent variable in S. Putting 
(fST)|) and (jSTj) in (|4*5| . one arrives at 

2V = 0(« a -i/ a ) + l, (52) 

which shows that (k — v) should be constant. One can use (|52|) in (l50l) and (IBTj) , 
to obtain 

a a ^a = -y [-^ («« + ^a) + 1], (53) 

&aC a = -y (« a + !]• (54) 



These two equations are not independent of (|47|) and (|48|) . The consistency 
condition is 

2 (Ka + ^) 2 -l=40 2 ?7 Q Q , (55) 

showing that [(« + v) 2 — A-qO] should be constant. Noting that (k — v) has to 
be a constant as well, the second condition can be stated as (k v — r\ 6) should 
be constant. So the conditions that the system be superautonomous are 

K a — v<x — constant, (56) 
K a v a — rj a 9 a — constant. (57) 

There are some special cases resembling those encountered in [7]. 

1 Constant coupling: Here H a does not depend on a, and A a tends to infinity 
(zero) if and only if the corresponding eigenvalue of Di is greater (smaller) 
than one. 
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2 Periodic coupling: 

H a +m = H a . (58) 

In this case the behavior of the eigenvalues of D a ^ is determined by the 
eigenvalues of D a + m ^ a : An eigenvalue of D a4l tends to infinity (zero) if 
and only if the corresponding eigenvalue of D a+rrita is greater (smaller) 
than one. 

3 Defects in the lattice: No new phenomena is seen, as long as the defects 
are localized, i.e. they are far from the boundaries. So if there is a lattice 
that has some defects but otherwise is uniform, the static behavior near 
the boundaries is similar to that of a uniform lattice |12j . 

4 A lattice with different behaviors at different end points: The behaviors 
of the static solution near the two ends are independent of each other, 
provided the behavior change occurs far from the boundaries. So all the 
phenomena seen in previous special cases can be seen at each boundary, 
independent of the other boundary. 

4 An example, the voter model 

The voter model is a lattice each site of which is either full (•) or empty (o). 
The reactions on a link are 

•o — y oo, with the rate u, 
•o — y ••, with the rate v, 
o» — > ••, with the rate it, 

o» — y oo, with the rate v. (59) 
Of course, the rates u and v may be link-dependent, using (|17l) . 

7]a V a , 
a 1i a , 

Vq. — Wa, 
(a =0, 

cr a = 0. (60) 



In order that the system be superautonomous, (1561) and (|5T[) must hold. (|56j) 

reads 

u a ~ v a = constant, (61) 
and (I5T1) is an identity. Assume that (fBTj) holds. Using (PHI) , one has 

1 

A := v . v . ■ (62) 
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From (|55|) . one has 

(j> 2 (u a - v a ) 2 = 1, 
one of the solutions to which is 

1 



Putting this in ([iTjl. ([15]). (|53| . and ((54)) , one arrives at 
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One set of solutions to these equations is 



u a 1 
v„ 1 



Putting this in (j44|) and (j45J) results in 

- v 



A,.= 



t — fl 







1 



One also has 



_ ,—1 



So 



One has 
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-1 l-4>v L 
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- 1 
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U a -1 
V a ■ 



tr(D aM ) = 1 



Vg-l ■ ■ ■ Vj-1 
U a -1 ' • ' U M +1 , 



showing that the eigenvalues of D a ^ are 

^ = 1, 

' «a-l ' • ' U/i+1 



A = 



Ua-1 • • -U M+ 1, 
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A special case is when 



(77) 



In this case the matrix E Q is singular. Yet one can obtain the matrix D afl as a 
limit (b — y oo of (1721). The result is 



1 



1 



V a -1 
1 



1 

Vl-l 



1 



1 



1 



1 



(78) 



and both eigenvalues of D afl become one. 
Let's study some special cases. 

1 Constant coupling: 

u a = constant. 

In this case A 2 tends to infinity if and only if v a is greater than u a . 

2 Periodic coupling: 



(79) 



(80) 



In this case the behavior of the eigenvalues of D a ^ is determined by the 
eigenvalues of D a + m>a , which are one and 



A 



2 ^/H-m— 1 ' " ' ^/x 

Ufi-^m — 1 ' " ' ^fj 



(81) 



If u a is greater than v a , then A 2 tends to zero. If u a is smaller than v a , 
then A 2 tends to infinity. 

3 closed lattice: Let the L'th site be the same as the O'th site. One has 

xl+^l = Xp. (82) 



Combining this with (1341) . it is seen that x M should be the eigenvector of 
Di,^„ „ corresponding to the eigenvalue one, which shows that 



(83) 



showing that the stationary profile of the density is uniform. This is 
despite the fact that the reaction rates are not necessarily uniform. 



5 Concluding remarks 

General autonomous exclusion models with nearest-neighbor interactions on a 
one-dimensional lattice were studied, for them the reaction rates were nonuni- 
form. By autonomous is is meant that the evolution equation for the expectation 
values of the number operators are closed. It was seen that the condition that 
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the system be autonomous is the same as the analogous condition for uniform 
lattices. A transfer matrix method was introduced to solve the equation for 
the static configuration of the expectation values of the number operators. The 
static phase picture of these systems, including possible phase transitions, was 
investigated. These phase transitions correspond to a reduction of the role of 
boundary reactions on the profile of the expectations of the number operators, 
similar to the case of uniform lattices. Also similar to the case of uniform 
lattices, these possible phase transitions are not affected by the boundary con- 
ditions. Moreover, all of the above mentioned possible static phase transitions 
near one boundary, are controlled by only the bulk reaction rates in a large part 
of the lattice one boundary of which is the same boundary. So phase transi- 
tions at different boundaries arc independent of each other, and finite defects 
far from boundaries have no effect on the phase transitions. It was seen that 
if the reaction rates satisfy ceratin additional conditions (which are essentially 
the constancy of ceratin combinations of reaction rates) a closed form can be 
obtained for the transfer matrix. Systems satisfying those conditions were called 
superautonomous. The example of the voter model was studied in more detail. 
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